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Abstract 

The general approach to chain equations derivation for the function generated by a 
Miura transformation analog is developing to account evolution (second Lax equation) 
and illustrated for Sturm-Liouville differential and difference operators. Polynomial 
differential operators case is investigated. Covariant sets of potentials are introduced 
by a periodic chain closure. The symmetry of the system of equation with respect to 
permutations of the potentials is used for the direct construction of solutions of the chain 
equations. A "time" evolution associated with some Lax pair is incorporated in the 
approach via closed t-chains. Both chains are combined in equations of a hydrodynamic 
type. The approach is next developed to general Zakharov-Shabat differential and 
difference equations, the example of 2x2 matrix case and NS equation is traced. 
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1 Introduction 

A potential reconstruction within some linear problem frame depends on a class in which one 
searches a solution of inverse problem [1]. The coefficients (potentials) of the equation itself 
naturally enter some algebraic structure generated by transformations that preserve the form 
of the equation operator. So the general variety of potentials splits to subsets invariant with 
respect to the action of the transformations. Such transformations as Darboux (Schlesinger, 
Moutard) ones (DT) are generated by transformations of eigenfunctions of a given differ- 
ential (in D) operator. The DT introduces naturally some intermediate object o = (D(f))<p 
that is linked to a potential by the generalized Miura transformation. The main ideas of 
G. Darboux, in principle, determine the form of the transformation [2] and the connection 
between potentials and a via a factorization of the operator under consideration [3, 4]. The 
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structure of the transformation depends on the ring to which the operator coefficients belong 
as well as abstract differentiation realization [5], [6]. The transformations and the natural 
Miura transformation for nonabelian entries (differential ring) are studied in [4] . Namely the 
object allows to operate effectively by the spectral problem data. We restrict ourselves here 
by one-dimensional problems, but important steps towards the multidimensions are made 
[7, 8]. 

If one substitutes the potentials in terms of the a for any iteration of the DT it yields in 
chain equations that are differential- difference equations and may be treated as the equiv- 
alent form of the spectral problem one starts with. It was discovered for the simplest one- 
dimensional problem in [9] and extended in [10]. The representation gives new possibilities 
to produce explicit solutions as well as study difficult questions of a potential uniform ap- 
proximation [11]. The technique is directly connected with the quantum inverse problem [12] 
and its use in soliton equation integration. In the next section we consider a general spectral 
problem, polynomial in a differentiation. We begin from the appropriate evolution equation 
and reduce the consideration to the stationary case that generates a spectral problem. Show- 
ing how the connection of the polynomial coefficients (potentials) with a appears, we derive 
examples of the dressing chain equation (DCE). 

The solution of a periodically closed DCE may be analyzed from the point of view of its 
bihamiltonian structure [13] as, for example, in [19]. Some important algebraic structures 
connected with the problem are studied in [21]. The symmetry of the system is naturally 
connected with the DT and generates a finite group that we use to simplify the problem. In 
the Sec. 3 we introduce projecting operators for the irreducible subspaces of the symmetry 
group and the corresponding variables, in the example we consider the transition to the new 
variables give a chance to express solutions of the DCEs in quadratures. 

The same function a satisfies an additional equation that appears if one studies the DT- 
covariance of the evolution equation of a Lax pair. A "time" evolution associated with some 
Lax pair is incorporated in the approach via closed t-chains, see also the Sec. 3. Both chains 
are combined in equations of a hydrodynamic type. 

The chain equations for matrix spectral problems as Zakharov-Shabat (ZS) one were 
considered in [14, 15], see also [11]. In the sections 5 we continue to study the general case 
and in 6 we treat an example of NS equation. The last section is devoted to the specific 
features of a DSE derivation for a difference ZS problem. 



2 The Miura maps and chain equations 

Let us consider the differential operator 

N 

L = Y,a n D n (1) 

n=0 

on a differential ring A, ip e A, with coefficients (potentials) a n e A and an evolution equation 

Ipy = 14 ■ (2) 

Here the operator D is a differentiation by some variable (or abstract one) and ip y is the 
derivative with respect to another one (see [6, 5] for details and generalizations). We also 
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would denote Dip = ip' . The transformation of the solutions of the equation is taken in the 
standard Darboux form 

iP[l] =DiP- aiP, (3) 

where 

a = 0V" 1 (4) 

with a linear independent invertible solution (p of (2). 

The transformation of the coefficients of the resulting operator 

N 

L[l] = a n [l]D n 

is defined by 

a N [l] = a N (5) 

and for all others n, by 

N 

a n [l] = a n + [akB k ,k-n + (a' k - aa k )B k ^ k ^ n ] (6) 

k=n+l 

that yields a covariance principle. It means that the function ^[1] is a solution of the equation 

iP y [l] = L[1M1]. 

The result is a compact reformulation of Matveev theorem [6]. 

The functions B„ hn are introduced in [4]. We reproduce here the definition and useful 
relations for these differential polynomials. 

Definition 

B n)0 (<7) = l, n = 0,1,2,..., 

and recurrence relations 

B n:k (a) = B n _ ljk (a) + DB^^a), k = 1, n - 1, n = 2, 3, . . . . 

B n>n (a) = DS n _i jn _i(<r) + B n {a), n = 1, 2, . . . . 
here and below the functions B n (a) are standard Bell polynomials. It yields 
k-i 



B n ,k(o-) = J2 [ n l k l +l ) B ™M) Dk ~ l ~ l °> k = l,n, n = 0,l,2,.... 



The formula gives the convenient algorithm for the evaluation of the generalized Bell poly- 
nomials. 

One results in the important conclusion: 
Statement 1. If a function a satisfies the equation 

D y a = o y = Dr + [r, a] (7) 
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where r = Y^n=o a n B n (cr); the operator L a = D — a intertwines the operators D y — L and 
D y -L[l\. 

Note that (7) the for a from (4) is the identity. 

For the derivation of the dressing chain equations we consider the stationary solution of 
the evolution equation (2) of the theorem: 

D y <j> = <M, (8) 

It gives D y o = 0. For example, in matrix case \i = diag{/ii, ...,//„}. Hence the corollary of 
(7) 

Dr + [r, a] = 0, (9) 

is the Riccati equation analogue, mentioned in the introduction, that we would name the 
generalized Miura map. It connects the "potentials" (coefficients of the operator L, see the 
expression for (1)) and a at every step (i) of DT iterations. Further we would supply the 
functions by the upper index i that show the number of iterations made. In the scalar case 
the commutator is zero and (9) reads r = c = const. If one have the only potential do 
with the rest invariant a l n = a^ +1 , n ^ 1, then the expression for the potential is a direct 
generalization to one mentioned in the introduction: 

N 

ai = -J2<Bn(vi) + d- (10) 

n=l 

In this case the derivation of the chain equation is made by the substitution of (10) into the 
last of (6) (n = 0), taken with the indices i,i+l at the correspondent side. The equation (10) 
in the case N=2 gives the link of the <jj with the potential Ui that enter the second order 
operator 

L = -D 2 + Ui . (11) 
Supplying the relation by indices, one have 

o-i, x + of + A*i = Ui- (12) 

Next nontrivial examples are connected with the case N=3 and the Miura transformation 
generalization (9) is given by 

3 

0-i, xx + -(0-i) x + of + UiGi +Wi = ^ (13) 

It connects the coefficients u^Wi with For both cases the DT has the same form 

u i+l = u i ~~ 2<7j jX , 

(for N=3 a 3 = 1, hence -2 should be changed to +3). Finally, if one starts from the second 
order Sturm-Liuville equation (11), its chain equation partner is 

o- i+ i, x + of +1 + = -o- i:X + o-] + (14) 
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For the third order operator if Wi = 0, the DCE is 
3 3 

(&i+l,xx + ^ ( a i+l)x + Cj+i - Hi+l)o~i = (<7i,xi + 2 (^i )x + ^ _ + "i® ix°~ i<3 i+\ (15) 

The case of zero «)i is obviously a reduction for the space of solutions of the linear problem 
and some modification of the DT formula for the eigen functions is necessary [23]. 

If a Lax pair of a nonlinear system is traced, one should consider one more operator, say, 

M 

A = J2KD n (16) 

n=0 

and the correspondent evolution 

A = Aip. (17) 

The coefficients of both equations (1,16) depend on a set of variables (potentials) u±, ...u n - 
eventual solutions of nonlinear equations. In the case of joint covariance property [5, 15] the 
DT transforms of the coefficients induce the DT of the potentials. The Statement 1 analogue 
accounts that the same a enter the operator L a 

Statement 2. If a function a satisfies the equations (7) and the equation 

D t a = a t = Dq + [q, a] (18) 

where q = J2^f=o b n B n (o~); the operator L a = D — a intertwines the pairs of operators 
D y — L, D t — A and D y — L[l],D t — A[l]. It means the integrability of the compatibility 
condition of the equations (2,17) in the sense of the symmetry existence with respect to DT 
of the potentials Ui — * 

The compatibility condition of the equations (7, 18) yields the extra equation 

Dq y + [q y , a] + [q, Dr + [r, a]] = Dr t + [r t , a] + [r, Dq + [q, a]] (19) 

That links the potentials and the element a. In the case of the unique potential it is possible 
to express the potential u as a function of a [16]. Considering the iterated potentials 

u % = fifa) (20) 

(now the index is again a number of iterations), allows to produce the dressing chain equation 
substituting the function into the DT formula Wj[l] = u i+1 . 

The scalar case is much more simple. From (19) it follows that D(q y — r t ) = 0, or 

M N 

%-n = J2 ( b ny R M) + K B' n {a)Dq) - £ (a nt B n (a) + a n B' n (a)Dr) = const. (21) 

n=0 n=0 

The good example of this case is the KP equation and its dressing chain, whence the potential 
is extracted from (21). In the theory of solitons (integrable equations) the case (15) generates 
Sawada-Kotera equation while (13) corresponds to the famous KdV. Other reductions are 
more complicated from the point of view of the chain equations derivation: it is need to 
express the potential from (13); e.g. so it is for the Boussinesq equation reduction case [23]. 
For the chain equation in this case see [16]. 
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3 Periodic closure and time evolution 



The periodic closure of the DSE (14) in the KdV case produces a finite system of equations 
that possesses the bi-Hamiltonian structure [13]. As the authors of the paper had written 
about the case N=3 (one-gap potentials) "It is a useful exercise to derive explicit formulas 
for <7j directly from equations of the chain ([9])." Below we briefly show how to do it and 
give the formula. There is an important question arises on this direct way: how to extract 
a potential dependence on the additional parameter t from a Lax pair for a compatibility 
condition? 

We propose to specify the "time" evolution via the t-dependence of x-conserved quantities. 
In this section we begin to study this problem in the terms of the same chain variables to be 
dependent already on the time. Let us start from the system for three functions <jj for the 
simplest nontrivial closure 

jz<7i{x) = cr 3 (x) 2 - a 2 (x) 2 + /i 3 - /x 2 , 

£a 2 (x) = a^xf- a 3 (x) 2 + m -na, (22) 
£<r 3 (x)= a 2 (xf - ax(x) 2 + /i 2 -/ii. 

The direct way in the bi-Hamiltonian formalism is initiated by the paper [13] via generating 
function and conservation quantities introduction. In this simplest example we have taken 
now the integrals solve the problem completely. If one express the third variable <r 3 as 
the linear combination of the rest ones by means of the first integral (Casimir function) 
c = <7i + <x 2 + <7 3 and substitute it into the rest equations of (22) then, after use of the second 
integral A, one arrives to the differential equation of the first order for elliptic functions. It 
is convenient to show this fact in the terms of other variables [13] gi(x) = cr 1 (x) + cr 2 (x), 
92{%) — <^ 2 (x) + cr 3 (x), g 3 (x) = <7 3 (x) + <J\{x). Let us exclude g 3 (x) by the connection 
9z{ x ) = 2c — g±(x) — g 2 (x). Further, omitting the argument x, the inverse transformations 
become 0\ = —g 2 + c, <j 2 = —c + gi +g 2 , <7 3 = c — gi Plugging the transforms into the system 
(15) we obtain two differential equations for the new variables 

4-01 = A*i - A*2 + 2c#i -g\- 2g 1 g 2 , ^-g 2 = -2cg 2 + 2g 1 g 2 + g\ + /i 2 - M3- (23) 
The second integral of motion in terms of gi is more compact 

^ = 919293 + ^203 + ^102 + A*301- (24) 

It allows to express g 2 as the function of g\ and one arrives at 



£gi{x) = -fa + 2cg x - g\ - 2(-g 2 + fi 2 - ^ + 2g ± c+ 
(gf - Agfc + 2gf(/i 1 + 2c 2 - 2/i 3 + ji 2 ) + 4(A - fj, 2 c - //ic)#i + fij + n\ - 2fi 2 fi 1 ) 1/2 



(25) 



The next problem is to solve the equation (25) by means of the theory of elliptic functions. 
The Weierstrass or Legendre canonical form of the integral yield a solution of the problem 
after Abel transformation [20] and use of the algebraic formulas that give <t 2)3 . Finally we 
have the explicit dependence of <7j on x and parameters c and A [17], see also Sec. 5. 

Let us turn to the problem of a time evolution arising from a Lax representation for 
some nonlinear equation. The main instruction is the search of the time dependence of the 
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x-independent entries c and A of the solution of the equation (25). Let us consider a DT- 
covariant evolution with the "time". In the KdV case the second Lax operator has the form 
(1). The account of a connection between the potential and a (12) produces the MKdV 
equation for the function a. The substitution into the equation the x-derivatives from the 
system (22) yields 

f (7 3 (X, t) = \ (/i! + /i 2 - 5/i 3 - F) ^CT 3 (x, t), . , 

F = erf + 2a 2 cr 3 + 1o x o 2 + o\ + 2a 1 a 3 + a\. y ' 

The equations for the a 2 .3 are written similar if one use the mentioned cyclic symmetry 
(in indices). Such equations are general; we would name them the t-chain equations. The 
form of the equations is typical for the so-called hydrodynamic-type equations . The system 
is diagonal and Hamiltonian, it could be integrated by the hodograph method [18]. The 
integration is trivial on a subspace of c=const, for F = c 2 . The final step of this direct 
construction is the use of the equation (26) and ones for <7 2 . 3 . First of all one can check that 

d 

—c = -3a?,fi 3 + 30-3/12 - 3o-3>i + 30-^3 - 3o-j> 2 + 3af^i = %Cx, 

that is zero if the x-dependence of the <7j is governed by (22). If one plugs the t-derivatives 
of (Tj from formulas like (21) into the derivative A t so that 

A == glgffin! - 3^3) + glgl(3/j,2 - 3//i) + g\g\{Znz - 3^ 2 ) + 

5 , 25 , i(6// 3 ^2-6//i//2)+5 , 25 , 3(6/i3^i-6//3^2)+5 , i5'3(6//i//2-6/i3Mi)+3 (/ii - /i 3 ) (a*2 - ^3) (a*2 - 

The analysis of the expressions leads to the statement 

Statement 2. If for all g i exists a common period X the x-independent polynomial c 
does not depend on t and A is linear function of t. 

For the proof it is enough to notice that the second order combination 5'2fi'i(6/i 3 /i2 — 
6/^1/^2) +fi , 2fi , 3(6/i 3 /ii — 6//3// 2 ) + gigs(6/iifj,2 — 6//3//1) is linear combination of x-derivative and 
constant, because, for example 

9i (.92 ~ 93) = ^92(x) + -^9s(x) + /i 2 - /ii 

After integration by x over the period X of the last equation and use of the conservation 
laws for the KdV and MKdV equations in combination with (15) one arrives to the linear 
dependence of AX on t. The coefficient may be recognized in the expression for A t : it 
is a combination of the eigenvalues /ij. For example, the / a^f af,i = 1,2,3 and other 
more complicated conserved quantities [25] should be accounted, the integrals over X of the 
combinations are zero. The resulting formula for A should be substituted into the solution 
of (25). 

Remark 1. It is shown (see (26)) that the derivatives of by t are linked with the 
derivatives by x so that the equations of hydrodynamic type appear. The equations could 
be integrated by the methods developed in [18]. 
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4 Discrete symmetry 



4.1 General remarks 

One can easily check the invariance of the every conserved quantity (e,g, mentioned here c 
and A) against the permutations of the elements a"j as well as the covariance of the systems 
(22) and (26). The symmetry with respect to the cyclic permutation of the variables (or 
indices) is obvious if one recalls its DT origin, hence, the observation is general. For the 
illustration we again exploit the simplest "KdV" case starting from the representation of 
[13, 21]. Let us consider the periodically closed DCE of the odd N=2g+1 and introduce a 
vectorial notations, namely 

S = (a 1 ,...,a N ) T , n = (//i, ...,n N ) T 

and 

L = (a x , ...,a N ) . 
The closed chain equation (14) may be rewritten either in the form 

or 

JV-l 

E*= £(-l) fe S fe (E 2 -/.). (27) 
k=i 

where the operator of permutation is represented as the matrix S [13]. Both forms are 
obviously invariant with respect to the S-transformation because the matrix S and operators 
in the equations (27) commute. 

The same statement is valid for the equations of the time evolution. 

Let us emphasize that the r.h.s. of the equations (27) is tensorial with respect to the 
components of S so the action of the group transformation is the tensor (direct) product of 
the group representations in correspondent vector spaces. 

If one introduce the cyclic permutation operators T s , its action determine the matrix S 

as 

T s CTi = Sik&k — &i+l(modN)i 

the powers of the matrix of the previous section produce the group, S k G C n C S n such that 
is valid for any symmetry group and give a basis for integration of the covariant equations. 
The technique uses the Poisson representation of the system (27) 

^^{Hft/j),^,^)}, (28) 

where the operator 

Hp,n) = j2(\<>i+Wd (29) 
k=i 6 

is invariant with respect to the group transformations and defines a linear operator adn with 
respect to the Poisson bracket 

{a l ,a k } = (-l) k ' i (l-5 tk ),k<t., (30) 
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It is easy to check, that 

n n 

{//(E), Oj} = ^(-KfT,j + /JfcW,ff 3 )) = 5Z(<Tfc + ^fc){ffc,0-j}, 
fc=l J fc=l 

that yields the system (27). 

The integration of the system may be understood in a "quantum mechanical" language, 
introducing first the "commuting" functions Cj(E) from the operator adn kernel Cj G X: 

a4C» = 0. (31) 

Next, the eigenvalue problem for a ^(E,//) outside the kernel can be considered as a matrix 
one at some basis 

ad H ipi = AjV'i, (32) 

i = 1, ■■■g,n — 2g + 1. The symmetry of the equations (31, 32) with respect to the transfor- 
mations T s (11) follows from obvious relation 

#(T S E,T S //)=#(E,//), (33) 

it means that matrices of a representation of the symmetry group commute with the matrix 
of adu on the correspondent subspace with constant Casimir operator. 



4.2 Irreducible subspaces 

The symmetry and the tensor structure of the r.h.s. of the equation (27) shows that the 
system may be simplified in the framework of Wigner-Eckart theory [26]. The statement of the 
Wigner theorem reads in the context of the equations (31,32) as the quasi-diagonal structure 
of the operator adu- Such structure and the correspondent simplification of the equations 
(31,32) solution is achieved by a transition to the basis of irreducible representations. The 
projecting operators pi to the irreducible subspaces are defined in the subspaces produced 
by chains that appears as the sum over group of the transformation action to some basic 
element. In the case of commutative group, chosen for a simplicity, the irreducible matrices 
are one-dimensional, and the basis is defined by the set of projectors 

Vl = Y J N l D\s)T s , (34) 

where iVj are normalizing constants, D l (s) are irreducible representations of the symmetry 
group and the T s is the group transformation operator in the space to be considered. 

In our case the operator coincide with the operator T s that have been introduced in the 
previous section and in the case of the cyclic permutation group C n the irreducible matrices 
are one-dimensional. Namely D j (e) = l,D j (s) = dj,D j (s 2 ) = a|, ... with a, = exp(j2m/N) 
where the integer N is the group order. Hence 

T sSj = a jSj . (35) 
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In the case of N=3 we project the system (23) or, originally (14) onto each subspace, 
having three equivalent equations. Let it be a = exp(2m/3); for example, the second of the 
resulting equations gives 

S2x = ni + an 2 + a 2 n 3 = 
°3 - °l + A*3 - fJ>2 + a (af - al + fxi - /i 3 ) (36) 
+a 2 {a\ - a\ + /i 2 - A*i) ■ 

Here denotes the r.h.s. of the equations (15). 

The inverse transform from original variables to the basis of irreducible representations 
reads 

o"i = s 1 + s 2 + s 3 , 
a 2 = si + as 2 + a 2 s 3 , 
cr 3 = s 1 + a 2 s 2 + as 3 . 



Plugging into the equation (36) yields 



d 



j^s 2 (x) = - (a - 1) ( 3a (s| + 2s!S 3 ) - afx 2 + - fi 2 + fi 3 ) , ^ 
£s 3 (x) = (a - 1) ( 3a (s| + 2sis 2 ) - a/i 3 + a/ii + /i 2 - /i 3 ) . 

The second conservation law (24) in terms of Sj allows to express the hamiltonian as a 
function of the only variable (s 3 ). The conservation laws are obviously the combinations of 
the irreducible polynomials o^. Together with the Hamiltonian (=A) conservation it leads to 
the spectral curve definition. 

Returning to general problem needs the tensor product space of the vectors S,/i. The 
problem of solution of the equations (31,32) simplifies if one use the mentioned symmetry, 
written in terms of 

N 

s t = Y, NtD^T^ = N- 1 ]T a\z\(J k . (38) 

sGG k=l 

Statement 3 By a direct application of the operator T s it could be checked that the tensor 
products of Si (see (35)) 

SiS^.-.Sj, (39) 

form a basis of irreducible tensors in the space of polynomials: the result differs from (39) by 
a constant factor aiCLk-.-Cbj. 

The further computations are conveniently made via the Poisson bracket 

{ 8j , Sl } = N- 2 ojliafrK **}■ (40) 

ik 

Particularly it is easy to show that the Cj could be chosen as a combination of irreducible 
polynomials (one could check this statement at the level of the examples from the mentioned 
paper [13]), the presented conservation laws are combinations of the irreducible polynomials 
of Oi, 
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Some adu- invariant chains are constructed by the following algorithm. We established 
that the first irreducible combination is proportional to the Casimir operator 

ad H si = 0, (41) 

Let us return again to the simplest example of N=3, taking the second equation of the basic 
system (e.g. (38)) 

ad H s 2 — n 2 — (a - 1) ( 3a (sj + 2sis 2 ) - a/i 3 + a/ii + // 2 - ^3) = (42) 
with the polynomial expression in the r.h.s.. Acting next, as 

adnn 2 = n' 2 = n^ 2 \ (43) 
one produces the higher order polynomial ar r.h.s; repeating, 

ad H n (l) = n (i+1 \ (44) 

one arrive at the statement which express 

00 

rf, = n w ,n (0) = s 1; = s 2 (45) 


as a polynomial function of s 2 and A. 

Statement 4 The eigenfunction ip of the operator ad H is expressed as the function of 
the variable A via solution of the algebraic (quadratic) equation 

= i) - s 2 - si. (46) 

The proof is delivered by the construction, taking into account the link (24) by which we 
express the powers s 2 with k < 3. 

An interesting link to the theory of automorphic functions applications to the system (27) 
could be found via the same symmetry. The example of (22) transformed as (37) allows to 
demonstrate it in almost obvious way. Let us recall that it is necessary to consider x as a 
complex variable. If one applies the operators T s or T s 2 to the both sides of the first and 
second equation of the system (37) and, simultaneously substitute 

x = a~ l x' , orx = a~ 2 x', (47) 

the equation rest invariant. In addition to the symmetry the double-periodicity of the solu- 
tions of the equations as elliptic functions (see also (32) or its Jacobi version by the transition 
■ipi = exp[Ci]) yields in Sj(x') = Si(x), 

x' — ax + r; a — 1, a, a 2 ; r = n\T\ + n 2 r 2 . (48) 

The property have place for any genus g, a e C 2g+ i. The solutions could be built in Poincare 
^-series [27]. See also [28]. 

Next, the t-dependence may be introduced via the scheme of Sec. 3. Other nonlinear 
systems are treated similar. Widen symmetry, for example from [21] includes reflections at 
the E, \i tensor product space and could give more information about solutions. 
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5 Oprerator Zakharov-Shabat problem 



Let us reformulate the general scheme of the dressing chain derivation in the nonabelian case 
from [15], starting from the evolution 

+ JDV = tf t , (49) 

with the polynomial first order L(D)-operator. The case is the nontrivial example of a general 
equation (2) with operator entries J and u (y is changed to t ). This way a form of evolution 
operator (Hamiltonian) is fixed as JD+u. Such form allows to consider one-dimensional Dirac 
equation [22] , or, in a stationary case, one goes down to a multilevel system interacting with a 
quantum field [24]. We would respect \1/ or other (necessary for DT construction) solution $ 
as operators. The equation could enter to the Lax pair of some integrable nonlinear equations 
as Nonlinear Schrodinger or Manakov ones. 

The potential u may be expressed in the terms of a from the equation (7), that we would 
rewrite as 

—a t + Jo x + [Jcr, a] = [a, u] (50) 

The structure of this equation determine the algebraic properties of the admissible dressing 
construction. 

The x-stationary version, when u x ,a x = 0,, the equation 50 yields 

-a t + [Ja,a] = [a,u] (51) 

that means (Spa) t = and the traceless possibility. The structure of a imply also the 
restriction 

Deta = DetM = JJ ^ (52) 
Namely, introducing the iteration index i, we have the link 

Du l + [u\ <Ji] = -DJ - [J, <7i]<7i, (53) 
the connection is linear, but contains the commutator. Let us denote ad a x = [cr,x\. Then 

«* = (ad ai - D)- 1 {DJ + [J, afa) . (54) 

The existence of the inverse operator in (54) need some restriction for the expression in () 
brackets, the expression should not belong to the kernel of the operator D — ad ai . In the 
subspace, where the Lie product is zero, the equation (53) simplifies. The DT is also simple. 

u l+1 = u i + [J, a t ] (55) 

Note that J is not changed under DT (see (5) Substituting the link (54) for i,i+l into (55) 
one arrives to the chain equations. One also could express matrix elements of u in terms of 
the elements of the matrix a and plug it into the Darboux transform (55) separately. 

Let us give more details of the construction in the stationary case, restricting DJ = 
and correspond to A, /i. Note that there are two possibilities for stationary equations 
from nonabelian (49): either 

^ t = \^ 
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or 

* t = (56) 

and the first of them leads to the essentially trivial connection between solutions and poten- 
tials from the point of view of DT theory [6] . 
In the second case one writes 

a = <Z> x <Z>- 1 = J- 1 ^®- 1 -u) (57) 

and DT takes the following form in terms of Sj 

tt^rVj-r 1 ^/], (58) 

where it is denoted s = here and further iteration number indices omitted. The 

potentials u 1 may be excluded from the equation (7) for this case 

a t = Dr+ [r, a], 

with 

r = J a + u = s. (59) 
The stationary case, after the plugging u l from (59) and returning indices, gives 

s i+1 = 8* + Ja i+1 - a\J (60) 

and 

Ds l +[s\a i ]=0. (61) 

links the derivative of s and the internal derivative of a. The formal transformation that 
leads to the chain equations is similar to the (54), namely, after substitution 

a 1 = -adJ 1 Ds i 

into the equation (60). 

Further progress in the development of this programm is connected with the choice of the 
additional algebraic structure over the field we consider. It can be useful for the concrete 
representation of solutions of the equation (61). For example, if the elements s l ,a l belong 
to a Lie algebra with structure constants C^p, then, after the choice of the basis e a one 
introduces the expansions (summation over the Greek indices is implied)) 

and 

Plugging into (61) gives the differential equation 

DC a + = o. 

If one defines the matrix 

B Pa = C^C V (62) 
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then, outside of the kernel of B 

rfr = -B£D& (63) 

By the definition of the Cartan subalgebra C the correspondent subspace does not contribute 
in the Lie product of (61). 

Statement. If, further, J belongs to a module over the Lie algebra, Je a = Jpa^p o,nd there 
exist an external involutive automorphism tau such that determine e a J (e.g. (ab) T = b T a T ). 
Then, the chain equation for the variables takes the form 

e +1 = & - B-^{D^)j, a + (jZpBrfDey, 

where the matrix B is defined by (62) and the components e a outside of C. Otherwise 

DC, = 0, 

if e 7 G C. Proof The statement is in fact the DT in the form of the equation (60) written 
in the basis of e a , in which the expression (63) is used. The subspace of C gives the second 
case. 

The system of differential equations is hence nonlinear as the matrix B depends on 
Remark. The scheme may be generalized for a nonstationary equation (2.5). The 
equation (5.2) then have the additional term Da\ from the r.h.s. 



6 Example of NS case 

The split NS equation is produced by the Lax pair based on (49) of the second order 

= (64) 

while the NS equations is solved under the reduction p = q. The choice J = <7 3 is obviously 
used. 
Let 

0~ = Viffi, u = v-iffi + U 2 C2 G sl(2, C) 
with the Pauli matrices <7j, as the basis e a 

[ai,a k ] = 2ie iks a s , (65) 

as generators of si (2, C). The matrix realization determine the left and right actions of J as 
usual matrix product, so the condition of the Statement of the previous section is satisfied. 
The "Miura" connection (51) is specified by 

Ja = ia 2 r] 1 - ia^ + r] 3 cr . (66) 

Plugging the result and u, a into the (51), one arrives at the equations: 

r][ + 2r} X r}z = 2ir) 3 u 2 , 

r]' 2 + 2r] 2 r] 3 = -2ir) 3 U!, (67) 
i] 3 — 2r)l — 2i]2 = 2ir] 2 ui — 2it)\U 2 . 
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From 52 it follows 

V3 = \t1M2 -Vi -V2, ( 68 ) 
that is in accordance with (67). The chain equations hence are 

u\ +1 = u\ - 2i4 , . 

4+ 1 = U 2 + 2iT][ { } 

where 

u i = -(V2/V3 + 2 ^2) /2i ( v 

«2 = (Vi/% + 2t7i)A 1 ; 

the correspondent iteration index % is implied and i] 3 is the function (68). The structure of 
the matrix elements of the potential u is such that q — U\ + i«2 so the reduction to NS case 
means the reality of both -Uj. The repulsive NS corresponds to p — —q*. The similar scheme 
is developed in [14]. 

6.1 The chain closures 

Simplest version of the closure, at the "first level' gives 
Introducing notations, 

vl = x, r]\ = y,vl = z, 
and taking into account the conditions (52) yield 



= x 2 + y 2 + z 2 , 
-fulfil = a\x 2 + a 2 y 2 + a\z 2 



providing restrictions for the constants 

11 2 2 / 22 22 2 2\/2 

-fi 1 fi 2 = x +y -{n^-a-tx -a 2 y)/a 3 , 



or 



, hence 

The equations (69) go to 
excluding z, one arrives at 



(72) 



2_ 2_ 2 1 1 2 22 

a 3 ~ a l ~ 0L 2i / i l/ i 2 Q; 3 ~~ / i l/ i 2 



CKl = ±tt3,Q;2 = ±«3- (73) 

(21 - l)x t /z = 2(a, + l)y, 

(f j -l)x t /z = 2(a l + l)x { > 



ln t (y) = ln t {x), 

1 + «2 l+a2 
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in a non-trivial conditions for constants (73), 

oti = -a 3 ,a 2 = -a 3 ,ln t (y) = -ln t (x), 

and, finally, gets 

y = c/x. (75) 
Plugging (75) into (74) and then z into (72) yields in 



\J — yUj/i^ 2 — X A — & 



= -(l-a 3 ) (76) 



that is again solved in elliptic functions [20]. The t-chains are obtained in a way similar to 
the section 3 using the second Miura map from Sec 2. 



7 Nonlocal operators 

Let again A be an operator ring, with the automorphism T. If for any two elements f,g<EA 

T(fg) = T(f)T(g), 

general formulas for DT for polynomial in T operators exist [6], see also [?] for examples. 
Here we continue to study the versions of ZS problem. We name the operator T as a shift 
operator, but it could be general as defined above. 

Let us take the general evolution equation in the case N=l. 

*p t (x,t) = (J+UT)^. (77) 

There are two types of DT in this case [?], denoted by indices ±. The DT of the first kind 
(+) leaves J unchanged. We rewrite the transform of U as 

U+ = a + (TU) (Ta+y 1 (78) 

where a + = 0(T0) _1 , further the superscript " + " is omitted. 

For the spectral problem correspondent to (77), the nontrivial transformations appear if 
in the stationary equation one introduces the constant element /i that does not commute 
with ip and a: 

(J + UT)cp = (pn. (79) 
The formula for the potential is then changed as 

U = wiTip)" 1 - Jo. (80) 

Let us derive the identity that links the potential U and a, doing it in a different manner 
than in [?] or, here, in the Sec. 6, starting from 

W 2 M = TM, (81) 
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and plugging it into the shifted equation (79): 



T(U)T 2 ( V ) = T( V )fi- JT(<p). (82) 

One has a Miura-like link 

aT(U)a = U +[J,a]. (83) 

T(a) = a^ 1 is accounted. Comparing with (78) yields new form of DT, that coincides with 
(84). 

U+[J,a]=U + . (84) 
Direct use of the equation (79) for expressing U in terms of r = (p[i(p~ l and a gives 

U = ra - Ja. (85) 

The element r is useful, also, for 

T{U)=a- l r- Ja' 1 . (86) 

Plugging (86), (85) into (83), one arrives at identity. The algorithm of the explicit derivation 
of the chain equations begins from the equation (83) solving with respect to U in appropriate 
way. For matrix rings, it may be a system of equations for matrix elements, that could be 
effective in low matrix dimensions of the "Miura" (83). 

The role of a + can play also the function s = ip^Tip) -1 . The equation (80) connects U 
and a. Let us rewrite (83) and the DT in terms of s, excluding U from (80) , denoting the 
number of iterations by index 

JJ[n\ = s n - Ja n . 

The equation (83) reduces to 

s n = <j n T{s n )(j n (87) 
The use of this result gives for the DT 

s n +i - s n = Jo- n+ i + a n J. (88) 

Then, solving the result (87) with respect to s one have the chain system. It could be made 
similar to the previous section by means Lie algebra representation. 

Let us mention that the chain equations for the classical ZS problem and two types of 
DT transformation were introduced in [?]. The closure of the chain equations specify classes 
of solutions. 



8 Conclusion 

Concluding, we wish to express a feeling that the technique elements we develop are general. 
Chain equation derivation is simply the result of substitution of a potential as the function of 
a into the DT formulas, but this problem of explicit form of the function could be non-trivial. 
The periodic closures of a chain for arbitrary N for KdV and other equations are studied very 
similar and leads to the expressions for the Oi and, consequently, for the potentials in hyper- 
elliptic functions by a construction. We also believe that the finite closures for the equations 
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may produce the solutions by the similar combination of a symmetry analysis for both x- and 
t- evolutions. The development of the technique for infinite chains do not look impossible as 
well. 

The work is supported by the Polish Ministry of Scientific Research and information 
Technology grant PBZ-Min-008/P03/2003. Author would like to thank M. Pavlov for fruitful 
discussion. 
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